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ABSTRACT 

Full sky surveys of peculiar velocity are arguably the best way to map the large scale 
structure (LSS) out to distances of a fewxlOOh^^Mpc. Using the largest and most accurate 
ever catalog of galaxy peculiar velocities Cosmicflows-2, the LSS has been reconstructed by 
means of the Wiener filter (WF) and constrained realizations (CRs) assuming as a Bayesian 
prior model the ACDM model with the WMAP inferred cosmological parameters. The present 
paper focuses on studying the bulk flow of the local flow field, defined as the mean velocity 
of top-hat spheres with radii ranging out to R - 500h 'Mpc. The estimated LSS, in general, 
and the bulk flow, in particular, are determined by the tension between the observational data 
and the assumed prior model. A prerequisite for such an analysis is the requirement that the 
estimated bulk flow is consistent with the prior model. Such a consistency is found here. 
At R = 50 (150) h^'Mpc the estimated bulk velocity is 250 + 21 (239 + 38) km s“'. The 
corresponding cosmic variance at these radii is 126 (60) km s ', which implies that these 
estimated bulk flows are dominated by the data and not by the assumed prior model. The 
estimated bulk velocity is dominated by the data out to R ^ 200h *Mpc, where the cosmic 
variance on the individual Supergalactic Cartesian components (of the r.m.s. values) exceeds 
the variance of the CRs by at least a factor of 2. The SGX and SGY components of the CMB 
dipole velocity are recovered by the WF velocity field down to a very few km s *. The SGZ 
component of the estimated velocity, the one that is most affected by the Zone of Avoidance, 
is off by 126km s '(an almost 2 sigma discrepancy). The bulk velocity analysis reported here 
is virtually unaffected by the Malmquist bias and very similar results are obtained for the data 
with and without the bias correction. 
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1 INTRODUCTION 


In the standard model of cosmology the large scale structure (LSS) 
of the universe grows out of a primordial perturbation field via 
gravitational instability. The continuity equation implies that the 
evolving density field is associated with a peculiar velocity field, 
both of which represents departures, or fluctuations, from a pure 
Hubble expansion. 

The Cosmic Microwave Background (CMB) dipole 
anisotropy, which is interpreted as the consequence of the 
peculiar motion of the Local group with respect to the CMB is 
the best evidence and example of that motion. First hints for the 
dipole anisotropy were given by |Conklm ^1969 | and |Hen 
and a more definitive determination by |Smoot et al.H1977k These 
discussions have set the stage for an effort that continues to the 
present epoch to map the three dimensional (3D) velocity field 
in our ’local universe’, namely our local cosmic neighbourhood 


(Tonry & Davis 

1981 |Aaronson et al. 1982| Lynden-Bell et al.| 

1988||Shaya et a 

l.|1995[ Willick et al.|1997[ Springob et al.|2007al 


|Tully et al.||2008| l. The motivation for these studies is twofold. 
One is the wish to map the underlying mass distribution in the 
local universe in general and the quest for the sources that induce 
the CMB dipole in particular. The other is to use the recovered 
local 3D flow field as a probe on the statistics of the underlying 
primordial perturbation field. This motivation leads immediately 
to the use of observed velocities as probes of the values of the 
cosmological parameters. 

The present epoch 3D velocity field is affected by non-linear 
processes, occurring mostly on small spatial scales. This has led 
people to focus on bulk velocities, namely the mean velocity of a 
large volume of space, as a mean of filtering out small scale non¬ 
linear motion and thus concentrate on the large scale modes that 
affect the bulk velocity (e.g. |Peebles|1980) l. The common choice is 
the definition of bulk velocity as the mean of a spherical top-hat 
window function convolved with the 3D velocity field. However, 
there are exception to this practice. The amplitude, scale depen¬ 
dence, rate of convergence and the direction of the local bulk ve- 
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locity are the subject of many studies (summary in |Watkins et al.| 

[ 200 ^ . 

Recent peculiar velocity studies have found contradictory re¬ 
sults in the amplitude of the cosmic bulk flow on scales of 50 - 
lOOh^'Mpc in excess of that expected in the A cold dark matter 
(ACDM) model ( [Watkins et al.|200^|Lavaux et al.|2010||Feldman| 
|et al.|2010| [Colin et al.|2011|( finding higher, seemingly inconsis¬ 
tent with ACDM, values while ( jNusser & Davis||2011[ |Tumbull| 
jet al.|20i2{|Ma & Pan|2014) finding lower values. Current surveys 
of peculiar velocities extend out to roughly lOOh^'Mpc, e.g. the 
SFI-(-(- survey ^Masters et al.|2006]|Springob et al.|2007b| l, and con¬ 
sequently current estimation of the bulk velocity do not go beyond 
that scale. Attempts to determine the bulk flow on larger scales are 
often based on using the anisotropic pattern of large-scale galaxy 
clustering, known as the redshift-space distortion ( |Kaiser|1987^ . A 
novel method based on the apparent dimming or brightening due to 
peculiar motion of galaxies in a redshift survey has been recently 
suggested jNusser et al.|2011) . 

Surveys of peculiar velocities are notorious for being ’diffi¬ 
cult data to use’. Observed radial velocities are noisy. The noise is 
driven by the uncertainties in the distance estimations. For exam¬ 
ple, the intrinsic scatter in the Tully-Fisher (TF) relation (jTu lly &| 
[Fisher] 1977| ( leads to a ~20% relative error in the distance, hence 
also in the estimated velocities. The relative error depends only 
weakly on the distance. Peculiar velocities datasets are sparse - the 
CF2 catalog consists of roughly 8000 data points within roughly a 
lOOh^'Mpc. The Galactic zone of avoidance (ZOA) limits the sky 
coverage of velocities surveys. The analysis of observational pecu¬ 
liar velocities datasets is hampered by various types of Malmquist 
bias. In the case of cosmic-flows2, the major one is due to the large 
errors on individual distance determinations. Such ’’difficult data to 
use” is best handled within the Bayesian framework of the Wiener 
filter (WF) and constrained realizations (CRs) of Gaussian fields 
( [Hoffman & Ribak|1991[[Zaroubi et al.|1995|[T999l l. 

We are using the largest ever and offered as the most accurate 
catalog of galaxy peculiar velocities, the Cosmicflows-2 dataset 
(CF2) ( [Tully et al.|2013] . The value of the Hubble constant used in 
the survey required for the decomposition of cosmic expansion and 
peculiar velocity components is calibrated by Type la supernovae 
calibration at z ~ 0.1 ( [Sorce et al.|2012] . The full 3D velocity field 
estimated from the CF2 dataset by means of the WF/CRs method¬ 
ology has been recently presented in [Tully et al.[j20i4t . The main 
aim of the present paper is to describe the flow field around us by 
means of its bulk velocity, namely the mean velocity obtained by 
convolving the full 3D velocity field with a given window function. 
Bulk velocities constitute a very compressed representation of the 
full 3D velocity fields and their analysis has played a leading role in 
the study of the cosmic flow field. We shall follow that tradition and 
dedicate the majority of the paper to the estimation of bulk velocity 
of the local universe. A secondary goal of the paper is the study 
of the motion of the Local Group (LG) with respect to the CMB. 
The problem to be addressed here is the assessment of how well the 
estimated LG velocity agrees with the observed CMB dipole. 

The structure of the paper is as follows. The CF2 data is de¬ 
scribed in ^ The WF/CRs methodology and its application to the 
CF2 data are reviewed in ^ ^presents a comprehensive review 
of the results, ^concludes the paper with a summary and a final 
discussion, 


2 DATA : COSMICFLOWS-2 

Cosmicflows-2 is the second generation catalog of galaxy distances 
built by the Cosmicflows collaboration. Published in [Tully et al.[ 
( [2013] , it contains more than 8,000 galaxy peculiar velocities. 
Distance measurements come mostly from the Tully-Fisher rela¬ 
tion and the Fundamental Plane methods jColless et al.[ j2001] ). 
Cepheids ([Freedman et al.[ ([2001f), Tip of the Red Giant Branch 
( Lee et al.H1993] ), Surface Brightness Fluctuation ([Tonry et al.[ 
( 2001] ), supernovae of type la ( [Jha et al!| ( [2007] ) and other miscel¬ 
laneous methods also contribute to this large dataset but to a minor 
extent (12%). 

For cosmic bulk flow purpose, in order to work above the scale 
of non-linearity in galaxy cosmological displacements (virial mo¬ 
tions in clusters, in knots of filaments, galaxy pairing interactions, 
etc..) we use in this paper the grouped version of cosmicflows-2 
containing 4,838 peculiar velocity values. Grouping also adds value 
to the following analysis since it decrease the statistical errors com¬ 
pared to individual galaxy values (~ 20%). The mean error for the 
groups is ~ 9% of the distance. See Tully et al. (2013) for details 
of the distribution of error estimates for components of the grouped 
catalog. 


3 METHOD 

3.1 Bayesian methodology 


The WF/CRs methodology has been extensively used to recon¬ 
struct the underlying density and three dimensional velocity fields 
from data bases of galaxy distances and peculiar velocities, starting 
with [Ganon & Ho ffman[ ( [T993] . A full Bayesian framework of the 
WF/CRs methodology was developed in [Zaroubi et al.[fl995] and 
its first application to peculiar velocities data base was presented 
by Zaroubi, Hoffman & Dekel (1999). The present WF/CRs recon¬ 
struction of the LSS follows the formalism of [Zaroubi et al.[ ( [l999] . 
The methodology was applied to the SEcat catalog of peculiar ve¬ 
locities jZaroubi[ ( 2002]) and m ore recently to the Cosmicflows-1 
data base ( [Courtois et ^.|2012 [[2013] and the Cosmicflows-2 data 
base jTully et al.|2014] . 

As mentioned earlier, the WF/CRs methodology recovers the 
linear density and 3D velocity fields. This is certainly a bad as¬ 
sumption for the density field but a good one for the velocity field in 
the quasi-linear regime, namely away from virial motions in groups 
and clusters of galaxies. As for the linear density field, the linear 
WF/CRs essentially recovers the divergence of the velocity field 
assuming the linear theory relation between the two fields. As the 
current paper focuses on the bulk velocity on scales of a few Mega¬ 
parsecs and above the linear theory provides the adequate theoreti¬ 
cal framework within which all present calculations are conducted. 


3.2 Prior model 

The standard model of cosmology, a flat A cold dark matter 
(ACDM) cosmology, is taken here as the prior model. The Wilkin¬ 
son Microwave Anisotropy Probe 5 (WMAP5) cosmological pa¬ 
rameters are assumed jKomatsu et al.[ ( [2009] ). The following pa¬ 
rameters, in particular, have been used here: Q,„ = 0.28 (the mass 
density parameter), h = 0.70 (the Hubble constant in units of 
100km s^'Mpc^') and erg = 0.817 (the r.m.s. of the linear den¬ 
sity fluctuations in a sphere of 8 h^'Mpe). The results reported here 
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remain virtually unchanged for the later release of cosmological pa¬ 
rameters of the WMAP team. The model further assumes that the 
primordial fluctuations constitute a Gaussian random field. 


3.3 Computational details 

Three kinds of velocities field are considered here: 1. A field ob¬ 
tained hy applying the WF to the Cosmicflows-2 data (V™); 2. A 
realization constrained hy the Cosmicflows-2 data 3. An un¬ 
constrained, i.e. random, realization The construction of a 

CR involves a one-to-one mapping from a random to a constrained 
realization and all the random realizations used here are the ones 
used to construct the CRs. The random realizations are constructed 
hy using the Fast Fourier Transform (FFT) technique. Hence they 
obey periodic boundary conditions. 

The different fields are constructed on a Supergalactic Carte¬ 
sian 256^ grid centered on the Milky Way galaxy. Three differ¬ 
ent computational boxes are used here with side length of L = 
320, 1280 and 2560h-‘Mpc (referred to as BOX320, BOX1280, 
etc.) The WF field has been constructed on the BOX320. The other 
larger boxes are used for the construction of the conshained and 
random fields. These large boxes are used to perform convergence 
tests on the bulk velocity of the CR fields. 

A note is due here on the role of the periodic boundary con¬ 
ditions and the construction of the WF, CR and RAN fields. The 
evaluation of the WF field at a given point in space relies on the 
calculation of various 2-point correlation functions. The analytical 
form of these functions is given by one dimensional (ID) semi¬ 
infinite integrals over the power spectrum convolved with functions 
such as the spherical Bessel functions ( [Zaroubi et al.|p^9^ . The 
integrals are evaluated by standard techniques of numerical ID in¬ 
tegrals. The value of the field at a given point is unaffected by the 
size of the computational box as the integrals used to calculate the 
correlation functions assumes an infinite space. The random real¬ 
izations, on the other hand, obey periodic boundary conditions and 
therefore are devoid of any possible power on scales larger than the 
scale of the box. As velocities are induced by long waves the size of 
the computational box reflects the tension between the needed res¬ 
olution (i.e. small scales) and the need to capture the power com¬ 
ing from long waves. The CRs are composed of the sum of the WF 
mean field and the random component and as such they are partially 
affected by the periodic boundary conditions. The large computa¬ 
tional boxes used here are chosen to minimize the effect and the 
convergence tests confirm it. 


3.4 Malmquist bias 


The study of issues associated with various biases that can arise in 
the measurement of distances to astronomical objects has a long 
history, dating from the pioneering studies of Gunnar Malmquist 


( |Malmquist|T920l 1 1 922| 1 1 924^ . 

In CF2 (Tully et al. 2013| l steps were taken to provide 
distances free of those that are often called the selection bias 
( |Willick||l994| ), but also confusingly referred to as simply the 
Malmquist bias. Even with statistically unbiased distances, how¬ 
ever, the forms of Malmquist bias due to scattering by errors from 
regions of greater sampling to regions of lesser sampling are ex¬ 
pected (whether of the so-called ‘homogeneous’ or ‘inhomoge¬ 
neous’ varieties). Also one must contend with the asymmetry in 
amplitude of peculiar velocity errors that arises from the translation 
of uncertainties that are symmetric in the logarithm of distances. In 


the following, it is these latter biases that will be addressed with 
reference to Malmquist biases. 

The WF reconstructed velocity field is used here to correct 
for the Malmquist bias. A full description and testing of the pro¬ 
cedure is to be published elsewhere and only a short description 
is presented here. The procedure is done in two steps. In the first 
the WF is applied to the raw data. New estimated distances are 
defined by subtracting from the observed radial velocity of a data 
point the radial component of the WF velocity field at its position. 
This adjustment brings the estimated distances to be much closer to 
their position in redshift space. The Wiener Filter is then applied to 
the revised catalog to yield the reconstructed flow that is presented 
here. The Wiener Filter correction of the Malmquist bias has been 
tested against mock catalogs drawn from N-body simulations con¬ 
strained by the CF2 data to reproduce the actual universe. 

The Malmquist bias introduces a spurious strong monopole 
term into the reconstructed velocity field but is expected to hardly 
alfect the bulk velocity which is associated with the dipole of the 
velocity field. Nevertheless all the analysis presented here is per¬ 
formed on the the Malmquist bias corrected data. A comparison of 
the bulk velocity determined from the corrected and uncorrected 
data is presented in Appendix [C| 

We have not attempted here to perform a maximum likelihood 
analysis to find the most likely ACDM model, but rather decided to 
choose an a priori WMAP5 set of cosmological parameters. This 
procedure leaves us with only one free parameter that enters the 
data covariance matrix, cr,, which serves as a crude proxy to small 
scales non-linear motion. We fine tune its value by the requirement 
that it makes the;^'“/d.o.f. equal to unity. Indeed for the CF2 data a 
cr, = 100 km s“' we find a,^^^/d.o.f. = 1.002 for the raw data and 
I.0I8 after applying the WF Malmquist bias correction. This small 
value of cr, is appealing as the data consists of a grouped version of 
the CF2 data base which largely negates non-linear motions associ¬ 
ated with groups and otherwise contains mostly spiral galaxies for 
which non-linear motions are quite insignifica nt. A comparison can 
be made with the cr, = 250km s“' adopted by Zaroubi et al. 1 2001 1 
in their study of a catalog of nearby early-type galaxies. In any case 
such a small cr, is relevant only for the determination of the very 
local velocities and it does not affect the bulk flow analysis. 


3.5 Bulk velocity determination 

Given the full three dimensional (3D) velocity field, evaluated on 
a Cartesian grid, the bulk velocity vector is readily calculated as 
the volume weighted mean velocity field within a top-hat sphere 
of radius R. The top-hat term is used here to denote a sphere with 
sharp cutoff. The bulk velocity is defined by: 

si I/"'"-’'*’' 

(Here, T stands for the type of the field - WF, CR or RAN.) 

Given that Vbuik(R) is obtained by a linear convolution over 
the full 3D velocity field, it follows that to the extent that v(r) is a 
Gaussian random field then Vbuik(R) is a Gaussian random vector. 
In the absence of constraining data, the prior expectation of of the 
bulk velocity is the null vector. For a large enough ensemble of 
random realizations the mean of Vbuik(R) over the ensemble is zero, 
(v*]'^'^*(R)) = 0. The mean field of the ensemble of CRs is the 

WF field, hence (v[,™'(R)^ = Given the finite number of 

realizations used here we expect some statistical deviations from 
the mean value. The variance of the dispersion of the ensemble of 
CRs and RANs from their mean values is used here as a proxy to 
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R [h-‘Mpc] 

WF 

CRs 

RANs 

10 

536 

543 ± 33 

511± 214 

40 

274 

282 ± 23 

322 ± 147 

50 

250 

258 ± 21 

292 ± 126 

100 

270 

279 ± 23 

209 ± 82 

150 

239 

241 ± 38 

161 ± 60 


Table 1. The bulk flow (in km s“*) of spheres of radius R of the Wiener Fil¬ 
ter reconstructed field and for an ensemble of 16 constrained (CRs) and ran¬ 
dom (RANs) realizations. For mean and standard deviations are presented 
for the two ensembles. 

the theoretical variance. In the following we use the term ’sigma’ as 
a shortcut for the standard deviation of the scatter around the mean. 
The statistical significance of the results will be measured by units 
of ’sigma’. 


4 RESULTS 

We start by presenting an overview of the WF reconstructed density 
and 3D velocity field. Figure[T]presents the WF density and veloc¬ 
ity field at the equatorial plane in supergalactic coordinates and the 
velocity field in the BOX320 in a 3D presentation. The Supergalac¬ 
tic Plane cut also contains the redshift space distribution of galaxies 
in an update of the 2MASS Extended Source Catalog ( [Jarrett et al.| 
|2000[ >. The 3D presentation of the flow field shows the dominant 
role of the Shapley supercluster in shaping the local flow. 


4.1 Bulk velocity 

Figure [^presents the main results of the paper. It shows the bulk 
velocity of the WF, 16 (BOX1280) constrained realizations and the 
matching unconstrained ones as a function of depth over the range 
of R = 10 to 150h“*Mpc. The compatibility of the estimated WF 
bulk velocity with the predictions of the ACDM/WMAP5 model is 
clearly manifested. The amplitude of WF bulk velocity lies within 
the expected scatter over the range of R = [10 - 150] h^'Mpc. At 
R = 150h“*Mpc the estimated amplitude is Vbuik = (239 ± 38) 
km s“*. It follows that we detect a bulk flow at more than a 6 
sigma significance at roughly the edge of the data. Some represen¬ 
tative values as a function of R are given in Table The lower 
panel of Figure shows the radial dependence of the three Su¬ 
pergalactic Cartesian components of the estimated WF, 

The plots also show also the Cartesian components of an ensem¬ 
ble of 16 (BOX160). At R ^ 150h-‘Mpc one finds: 

Vtaik [SGX, SGY, SGZ] = [-196 ± 43,59 ± 27,-123 ± 41]km s^'The 
flow is mostly within the Zone of Avoidance (ZOA) which almost 
coincides with the SGX-SGZ plane at SGY=0. The smallest com¬ 
ponent of the bulk flow is in the SGY direction. The largest com¬ 
ponent is in the negative SGX direction, detected with a 4.5 sigma 
significance. 

Using the largest computational box, BOX2560, the mean and 
the scatter around the mean of the amplitude of the bulk velocity of 
an ensemble of 20 CRs and their corresponding RANs are shown 
in Figure]^ all the way up to R = 600h“*Mpc. Some representa¬ 
tive values are presented in Table Again, the mean bulk veloc¬ 
ity lies within the 1 sigma bracket of the assumed prior model. At 
R = 600h-‘Mpc the WF yields Vbuik = (59 ± 18) km s-', a 3.2 
sigma detection. The expected amplitude for random realizations is 


R[h-‘Mpc] 

WF 

CRs 

RANs 

50 

11 

9.6 ± 6.6 

39 ±39 

100 

7.5 

6.6 ±5.4 

53 ±40 

150 

16 

13 ±9.2 

59 ±36 


Table 3. The alignment (in degrees) of V|,^Jj,(R) with V^^J]^(10h“*Mpc). 
Second column present the alignment for the WF field. The third and fourth 
columns show the mean and standard deviation for the ensemble of 20 con¬ 
strained and random realizations. 

(52 ± 28) km s“*. The comparison of the individual Cartesian com¬ 
ponents is even more illuminating. Considering the limited number 
of realizations employed here we expect to find statistical scatter in 
the values of mean quantities. For random realizations one expects 
the mean of the Cartesian components of the bulk velocity to be 
zero and the scatter to be isotropic. The deviation from these ex¬ 
pectations enables one to estimate the statistical uncertainty in our 
results. At R = 600h“'Mpc the mean over the ensemble of RANs 
of the individual components is a very few km s“'. We find here a 
non vanishing, at the 1 sigma level, SGX/SGY/SGZ component of 
the bulk velocity out to R x 400/200/300 h^'Mpc. 

4.2 Alignment 

The direction of Vbuik remains coherent over a large range of scales. 
Is it consistent with the ACDM/WMAP5 model? This question is 
addressed by calculating the angle between Vbuik at a given ra¬ 
dius and at R = lOh^'Mpc. Figure shows the angle between 
'^buik(10 h-‘Mpc) and V™^(R) for T= WF, CR and RAN as a func¬ 
tion of R. Table shows some representative numbers for the WF 
field and the mean and scatter for an ensemble of 20 constrained 
and random realizations. Tableand Figure]^ show that the bulk 
velocity vector can zigzag quite drastically in unconstrained real¬ 
ization but nevertheless some finite fraction of the realizations re¬ 
tain their coherent direction, akin to our own local flow field. 

4.3 CMB dipole 

In the standard model of cosmology the CMB dipole anisotropy 
is induced by the Local Group velocity with respect to the CMB 
background. This velocity should be recovered as the limiting case 
of the bulk velocity of a top-hat sphere of a null radius. This is ham¬ 
pered by two possible effects. One is the incomplete sky coverage 
of the velocities induced by the Zone of Avoidance and the other is 
possible non-linear effects that might affect the motion of the Local 
group. Table [^presents the LG velocity, as inferred from the CMB 
dipole anisotropy ( |Tully et al.|2008[ and references therein) 

The comparison shows that the WF recovers the SGX and 
SGY components of the CMB dipole velocity remarkably well and 
somewhat less well, at L7cr, away from the SGZ component. The 
deficiency of the WF reconstruction in the SGZ direction is ex¬ 
pected - it is the direction that lies deep in the ZOA and closest to 
the Galactic center. The surprising result is how well the WF flow 
field recovers the other two components, and in particular the SGX 
component that lies in the Galactic plane. 

The misalignment of the LG velocity of an ensemble of 15 
CRs with the observed dipole is found here to be (8 ± 7)°. However, 
limiting the study of the alignment in the SGX-SGY plane, namely 
the Supergalactic Plane, and thereby excluding the SGZ component 
a much better alignment is found, of (0.4 ± 5)°. 
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Figure 1. Left panel: The WF 3D velocity field (represented by vector arrows) and the linear density field on the Supergalactic Plane. Red contour denotes 
the zero over-density and the black contours the positive over-densities. Galaxies from the 2MASS redshift survey are presented by green dots and are placed 
in redshift space. Right panel: The WF of the 3D velocity field (represented by flow lines) in a box of L = 320h“*Mpc centered on the LG. The density field 
on the Supergalactic Plane is presented as well. 




CRs 




RANs 



R 

Fbulk 

Vx 

V'y 

Vz 

^bulk 

Vx 


Fz 

40 

309 ± 29 

-191 ± 30 

229 ± 35 

-64 ± 34 

318 ± 142 

-12 ± 178 

19 ±212 

40 ± 218 

120 

276 ± 30 

-195 ± 33 

123 ± 29 

-144 ±40 

254 ± 113 

-6 ± 127 

30 ± 138 

20 ± 116 

200 

201 ± 45 

-167 ±54 

48 ± 35 

-83 ± 37 

150 ±61 

-9 ±97 

17 ± 103 

17 ± 82 

320 

115 ± 55 

-91 ± 62 

19 ±46 

-30 ± 32 

101 ± 42 

-1 ± 71 

9 ±70 

15 ±47 

400 

90± 37 

-60 ± 55 

12 ±42 

-19±31 

80 ±34 

4 ±61 

6 ±53 

10 ±34 

520 

70± 22 

-36 ± 45 

6±31 

-14± 30 

61 ± 31 

5 ±50 

2 ±38 

5 ± 28 

600 

59 ± 18 

-27 ± 39 

5 ± 26 

-11 ± 29 

52 ±28 

6 ±43 

2±31 

4± 27 


Table 2. The mean and standard deviation of the amplitude and the SGX, SGY and SGZ components of 20 CRs (columns 2-5) and their corresponding RANs 
(columns 6-9). (CRs and RANs are calculated within BOX2560.) 


4.4 Statistical significance 


The standard deviation of the Cartesian components of the bulk 
velocity of the the ensemble of 20 constrained and random real¬ 
izations is shown in Figure]^ BOX2560 is used here to calculate 
the scatter out to R = 600h“'Mpc. The plot shows that out to 
R = 400h“*Mpc the variance in the value of each of the Cartesian 
components of the CRs is smaller than cosmic variance. The con¬ 


strained variance converges to the cosmic one on larger scales, with 
the rate of convergence varying with the Supergalactic directions. 

4.5 Comparison with SFI-n- and the COMPOSITE datasets 

[Nusser & Davis| ( [2011^ estimated the bulk flow of the local ve¬ 
locity field from the SFI-l-l- survey. They used ACSE model (see 
below), which is close in its approach to the WF/CRs, and their 
definition of the bulk flow is identical to the one used here. A de- 
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Figure 3. Upper panel: Amplitude of the bulk velocity: The mean (thick lines) and the mean ± one standard deviation (thin lines) of an ensemble of 20 
constrained (solid lines) and unconstrained (dashed lines) realizations. Lower panel: Mean (thick lines) and the mean ± one standard deviation (thin lines) of 
the three Cartesian components of the bulk velocity taken over the ensemble of 20 CRs (BOX2560) over the range R = [0 - 500] h^'Mpc: SGX (full lines), 
SGY (dotted lines) and SGZ (dashed lines) 


tailed comparison of |Nusser & Davis| ( |2011| l results is shown in 
Table which shows a close agreement. The CF2 reconstruction 
has a tighter error bars, by almost a factor of 2, compared with the 
|Nusser & Davis| ( |201 1| > numbers. A discussion of the comparison 
is presented below, the differences are mainly due to the improve¬ 
ments in the observational data. CF2 is more extended, consists 
of a larger sample and has a typical statistical fractional error of 


~ 18% compared with the ~ 23% of the SFI-H-. statistical error of 
~ 23% of distance on individual data-points, while Cosmicflows- 
2 grouped version (which includes some of the best observational 
raw data available in SFI+-t-) is somewhat smaller (~ 18%). 

The COMPOSITE compilation ( [Feldman et al.|201^ is anal¬ 
ysed here in the same way as the CF2 data. Namely, the bulk mo¬ 
tion of a top-hat sphere of the WF and CRs flow field. Resulfs for 
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R[h-‘Mpc] 

Data 

Vbuik 

Galactic (/, b) 

Vbuik,X 

^bulk.y 

Vbuik ,2 

40 

CE2 

274 ± 23 

oo 

00 

CM 

-226 ± 23 

102 ± 20 

-116 ±28 


SFI-l-t- 

333 ± 38 

(276°, 14°) 

-242 

103 

-203 

50 

CE2 

250 ± 21 

(280°, 18°) 

-209 ± 22 

90 ± 19 

-100 ± 24 

COMPOSITE 

243 ± 28 

(284°, 18°) 

-193 

89 

-1183 

100 

CE2 

270 ± 23 

(283°, 19°) 

-211 ± 31 

105 ± 23 

-131 ±33 


SFI-t-t- 

257 ± 44 

(279°, 10°) 

-198 

61 

-151 

COMPOSITE 

300 ± 33 

(282°,11°) 

-240 

76 

-163 


Table 5. A comparison of the bulk velocity obtained from the WF/CR reconstruction of the CF2 data and of the COMPOSITE data compilation and the 
jNusser & Davis|201H analysis of the SFI++ data. The bulk velocity of a top-hat sphere of radius R is presented in galactic (/, b) and in Supergalactic 
Cartesian coordinates. Velocities are in units of km s“*. 



VcMB 

VcMB.x 

VcMB.y 

VcMB.z 

Observed 

631 ±20 

-381 

331 

-380 

WF/CRs 

554 ± 56 

-364 ± 56 

321 ± 50 

-264 ± 70 


Table 4. The CMB dipole velocity (norm and the SGX, SGY and SGZ 
components all in km s“*): A comparison of the velocity evaluated at the 
position of the LG with the observed one. The WF/CRs denotes the esti¬ 
mated velocity of the Local Group. The one sigma scatter is the standard 
deviation taken over 20 CRs. 

R = 50 and lOOh^’Mpc are presented in TableAgain, a good 
agreement is found with a tighter error bars for the CF2. 

4.6 Comparison with Cosmicflows-1 dataset 

Cosmicflows-2 is the second generation catalog of galaxy distances 
built by the Cosmicflows collaboration. It has been preceded by the 
Cosmicflows-1 (CFl |Tully et al.|2008) database, which consists of 
data confined predominantly to redshifts smaller than 3000km s“*. 
It follows that the CFl is a subset of the more extended CF2 dataset. 
The |Courtois et aL] ( |2012| l study of the CFl database was a precur¬ 
sor of the present paper, using essentially the same methodology at 
the one used here. Given the shallowness of the CFl data lCourtoisI 
|et al.| ( |20T^ f ocused their efforts mostly on the local cosmography. 
Appendix |D| presents a comparison of the bulk velocity estimated 
from the CFl and CF2 data bases. The comparison serves as a man¬ 
ifestation for the power of the WF/CRs Bayesian methodology and 
it demonstrates how to interpret the results of the WF/CRs analysis. 
The comparison shows the following expected facts: i. the uncer¬ 
tainty in the estimated CF2 bulk velocity is much smaller that those 
of the corresponding CFl one; ii. the CF2 estimated bulk velocity 
is roughly within the one standard deviation scatter around the CFl 
estimated bulk velocity. 


5 SUMMARY AND DISCUSSION 

The paper reports on the application of the Wiener filter and con¬ 
strained realizations methodology to the Cosmicfiows-2 data set 
and presents an estimation of the cosmological bulk flow. A bulk 
flow is defined here as the mean velocity of a top-hat sphere cen¬ 
tered on the Milky Way. The bulk velocity has been studied by 
means of its amplitude, its three Supergalactic cartesian compo¬ 
nents and the alignment of its direction across different distances. 
In the limit of vanishing top-hat radius the bulk velocity should con¬ 
verge to the CMB dipole. This condition is indeed the case, with a 
caveat. The reconstructed SGX and SGY components of the LG 


velocity essentially coincide with the observed values. The SGZ 
component which lies deep in the ZOA and is most affected by it is 
roughly two sigma away from the observed value. On a larger scale, 
of R = 100 h^'Mpc, the amplitude of the bulk flow is 270 ± 23 km 
s“*. A non-vanishing signal is discovered all the way to the edge of 
the reconstruction with Vbuik = 59 ± 18 km s“' at R = 600h“‘Mpc. 
The behaviour of the bulk velocity, in terms of direction and am¬ 
plitude, is consistent with the predictions of the ACDM/WMAP 
cosmological model on all scales, ranging from zero radius out to 
R = bOOh^'Mpc. The consistency of the estimated bulk flow with 
the predictions of the assumed Bayesian prior model constitutes a 
prerequisite for the present analysis. Indeed, such a consistency has 
been found here. 

The Bayesian framework employed here is standing on two 
pillars. One is the WF which provides the mean, and also the most 
probable, field given the data and the prior model. The other is the 
ensemble of CRs, which samples the variance of the scatter around 
the mean field. There is tension between the data and prior model. 
The WF represents the equilibrium between these two components 
- the WF solution is dominated by the data where the data is strong 
and by the prior model where the data is weak. The CRs sample 
the variance around the WF. The stronger is the data the smaller is 
the variance of the scatter around the WF field, compared with the 
unconstrained cosmic variance. It follows that a comparison of the 
variance of the CRs with the cosmic variance provides a measure 
of the robustness of the reconstruction. The test is best manifested 
in Figure]^ which presents the variance of the three Supergalactic 
Cartesian components of the bulk velocity as a function of depth. 
As expected, the SGZ component is the one that is least constrained 
by the data and hence the constrained variance of this component 
converges to the cosmic variance at 400h“*Mpc. For the other 
two components the inequality extend to beyond R = 600h“‘Mpc. 
On the scale of R = 200h“'Mpc the estimated bulk flow is clearly 
dominated by the data, as is manifested by the fact that for Carte¬ 
sian components of the bulk flow the cosmic variance is larger than 
the constrained one by at least a factor of 2. 

Two recent studies have focused on the evaluation of the bulk 
flow of the local velocity field. [Nusser & Davis| ( |2011| l analysed 
the SFI-t-t- data base by means of the All Space Constrained Esti¬ 
mate (ACSE) model. The model seeks to reconstruct the full space 
3D velocity field using a Bayesian regularised fitting of base func¬ 
tion sampled from a prior model. Given such a full space recon¬ 
struction the bulk flow is determined as the spatial average over a 
spherical top-hat window. The ACSE model is akin to the current 
WF/CR method in its use of a Bayesian prior model and definition 
of the bulk flow, hence its results can be directly compared with 
the present ones. Table presents a comparison of the bulk flow 
determined from SFI-n- by the ACSE method and from the CE2 
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WF/CRs analysis. There is a good agreement between the two es¬ 
timates to within one or two standard deviations of each method. 
Given the different methods of analysis and the different data used 
such a discrepancy is expected. In both cases the estimated bulk 
velocity is compatible with the cosmic variance predicted by the 
WMAP ACDM model. 

[Feldman et al.| ( [20T0t used a different approach to the problem. 
They devised an ’optimal minimum variance’ weighting scheme 
and applied it to the COMPOSITE catalog which consists of a 
compilation of different peculiar velocities data bases. The Feld¬ 
man et al method does not provided an estimate of the bulk ve¬ 
locity of a top-hat sphere of a given radius but rather a bulk ve¬ 
locity on an effective scale of an optimally designed window func¬ 
tion. It follows that the Feldman et al estimated bulk velocity can¬ 
not be directly compared with the ones reported here and with 
[Nusser & Davis] j2011^ . [Feldman et al.| p010[ l reports a bulk ve¬ 
locity on scales of ~ lOOh^'Mpc of a magnitude of 416 ± 78 km 
s^'towards (galactic) {l,b) = (282° ± 11°, 6° ±6°) or (supergalactic) 
(sgl, sgb) = (168°, -34°). It is further reported that “[t]his result is 
in disagreement with cold dark matter with ... WMAP5 cosmolog¬ 
ical parameters at a high confidence level”. The incompatibility of 
the Feldman et al estimated bulk flow with the WMAP5 ACDM 
model stands in a sharp disagreement with [Nusser & DavisH2011) 
and with our results. This incompatibility is due to the different 
mode of analysis of [Feldman et al.[ ( |2010[ l and not due to their data 
compilation. The application of the WF/CRs methodology to the 
COMPOSITE dataset results with a bulk flow in full agreement 
with CF2 and with the |Nusser & Davis] ( [2011^ analysed the SFI-l--(-. 
In particular the WF/CRs bulk flow of the COMPOSITE dataset at 
R = lOOh^'Mpc is Vbuik = 300 ± 33 in the direction of (galactic) 
(/, b) = (282,11)° or (supergalactic) (sgl, sgb) = (163°, -33°). 

The standard interpretation of the observed dipole moment of 
the CMB temperature anisotropies is that it is the outcome of our 
motion with respect to the CMB. After making a suitable trans¬ 
formation to the LG frame of reference the CMB dipole is being 
expressed as the LG motion with respect to the CMB, denoted here 
by VcMB- The desire to understand the sources of that motion and 
how it integrates into the large scale velocity field is one of the main 
motivation for studying the large scale flow field. The current study 
share this quest. An on-going debate in the community is on the role 
of non-linear dynamics in shaping the CMB dipole (e.g. |Chodor-| 
[owski|201 l[|Nusser et al.|2014[ ). Our analysis sheds some light on 
the problem. The WF/CRs reconstruction is done within the linear 
theory framework. Any possible discrepancy between the LG ve¬ 
locity of the WF/CRs flow and the observed one can be attributed 
to limitations of the method, shortcomings of the data and non¬ 
linear effect. The present WF/CRs analysis reconstructs the SGX 
and SGY components of the CMB dipole remarkably well. This is 
not the case with the SGZ component, with an almost 2 sigma dis¬ 
crepancy between the prediction and observation. The SGZ compo¬ 
nent lies within the ZOA close to the direction of the Galactic center 
and thereby it is the Supergalactic component that is most affected 
by the ZOA. The good agreement in the SGX and SGY directions 
and the discrepancy in the direction of the Galactic center suggest 
a minimal role of non-linear dynamics and that the discrepancy be¬ 
tween the linear theory estimation and observation is dominated by 
the obscuration within the ZOA. 

The bulk velocity statistics constitutes a powerful statistics 
and cosmological probe of the large scale structure. The Bayesian 
methodology of the Wiener filter and constrained realizations has 
been used here to provide a robust estimation of the bulk veloc¬ 
ity of top-hat spheres on scales ranging from the CMB dipole up 


to roughly 200h“‘Mpc. The reconstruction has recovered the ob¬ 
served CMB dipole, away from the direction of the Galactic cen¬ 
ter, suggesting that the CMB has been hardly affected by non¬ 
linear processes. The resulting bulk velocity is consistent with 
the Bayesian prior ACDM/WMAP5 standard model up to roughly 
200h-'Mpc. 
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APPENDIX A: WE, CRS, AND ALL THAT 

The derivation here presents the essence of the WF/CRs reconstruc¬ 
tion. The reader is referred to [Zaroubi et al.[ ( [T999[ l for a full deriva¬ 
tion and for explicit calculation of the covariance matrices. In the 
following an ensemble of observed peculiar velocities is assumed, 
ju/j ^ ^ is assumed. It is further assumed that the radial compo¬ 
nent of the velocity is measured at its true distance, r,, i.e. the se¬ 
lection Malmquist bias is nulled. An individual data point is related 
to the underlying density field by: 

u° = v(r,) • f, -I- e, = M,- + Si. (Al) 

Measurement errors are assumed to be statistically uncorre¬ 
lated and the errors covariance matrix diagonal and is written as: 

(o'? + cffjSij (A2) 

Here cr, is the direct measurement error that typically scales with 
the distance of the data point. A cr, term is introduced to the er¬ 
ror covariance matrix as a crude proxy for the contribution of non¬ 
linear small scale motions. 

The estimated WF velocity field is obtained by a linear convo¬ 
lution on the data vector, |m”|: 

v'^'^fr) = (v(r)u°)(t<”M”) (A3) 

Given a random realization of the underlying velocity field 
a mock data base is constructed , i?,. This is used to construct a 
constrained realization: 

vCR(r) = v(r) -t ly(r)uf)(u‘’u‘’)^\u) - u-,) (A4) 

Given a prior model for the underlying velocity field and an 

errors model, the data auto-covariance matrix is written as: 

Rij = (m>”) = + (e/ey) = f,(v(r,)v(ry))fy -i- (err + crj)(5,y. 

(A5) 

The reduced of the data given a prior model is defined by: 


The value of cr, is tuned to yield = 1.0. 
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APPENDIX B: CONVERGENCE TEST 

Random realizations of the velocity field are used to construct con¬ 
strained ones. With the [Hoffman & Ribak| ( | 199 ft method there is a 
one to one relation between the random realization (RAN) and the 
resulting constrained one (CR). RANs are generated by means of 
FFT, which imposes periodic boundary conditions. It follows that 
the RANs do not have any power from outside the computational 
box. Two boxes are used here, BOX1280 and BOX2560. Figure 
|B 1 1 presents the amplitude of the bulk velocity of an ensemble 20 
CRs constructed within these two boxes. Note that the realizations 
used in each box are different. The distribution of both boxes are 
virtually the same and differences are consistent with the scatter 
expected for an ensemble of 20 realizations. 


APPENDIX C: BULK VELOCITY AND THE MB 
CORRECTION 

Table [CT] presents the bulk velocity, Vbuik(R), for the WF veloc¬ 
ity field calculated from the raw data and the Malmquist cor¬ 
rected data, evaluated at some representative depths. The compari¬ 
son shows a very small difference in amplitude at R = lOh^'Mpc. 
For R = 50h“*Mpc and larger radii the differences in amplitude 
amount to a very few km s“*, and the differences in the direction 
are of the order of a very few degrees. The differences between the 
outcome of the raw and Malmquist bias corrected data are much 
smaller than the expected scatter around the WF mean field. 


APPENDIX D: BULK VELOCITY: COMPARISON OP 
THE CEI AND CE2 DATABASES 

The present appendix aims at demonstrating the nature and power 
of the WF/CRs methodology by considering the estimation of the 
bulk flow from the Cosmicflows-2 (CF2) and the Cosmicflows-1 
(CFl) databases. The WF provides a Bayesian estimation of an un¬ 
derlying (velocity) field given the data, its measurement errors and 
an assumed prior model. The resulting estimation depends on the 
tension between the quality and strength of the data and the as¬ 
sumed prior model. Where the data is strong the recovered field is 
expected to follow closely the data, otherwise the estimation would 
tend to the null hypothesis of the prior model. The tell-tale sign of 
the data dominance is the variance of the constrained realizations 
around the WF field, which we call here the constrained variance. 
The constrained variance For a data dominated constrained field the 
constrained variance is smaller than the cosmic variance, namely 
the variance of unconstrained realizations. It is the constrained vari¬ 
ance which gauges the role of the data and not the value of the WF 
estimated field. 

CF2 contains as a subset the CFl database, which extends 
roughly out to 30h“*Mpc. The CF2 is a better database in the 
sense of its depth and number of data points. The comparison of 
the WF/CRs reconstruction out of the the CFl and CF2 datasets 
demonstrates how the quality of the data affects the constrained 
variance. Figure [DT] shows the WF (mean) bulk velocity as a func¬ 
tion of depth, Vbuik, and the ± scatter around the mean of the CRs 
for the CFl (in blue) and CF2 (in black) datasets. The case of un¬ 
constrained realizations is shown as well (in red). Two main points 
emerge from the comparison: a. Indeed the CF2 constrained vari¬ 
ance is considerably smaller than the CFl constrained variance and 
of the cosmic variance (of the random realizations). The CFl looses 


its constraining power, in the sense of it constrained variance con¬ 
verging to the cosmic variance, at about R ~ SOh^'Mpc. For the 
CF2 this occurs at R ~ 300h“*Mpc (Tablej^. b. There is no incon¬ 
sistency between the WF/CRs reconstruction out of the CFl and 
CF2 data. The CF2 mean (WF) Vbuik lies over all scales within the 
scatter of the CFl estimated bulk velocity. 
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R [Mpc/h] R [Mpc/h] 


Figure Bl. A convergence test on the size of the computational box: Plots show the bulk velocity calculated from an ensemble of 20 CRs calculated within 
the BOX1280 and BOX2560. Left panel: The bulk velocity of the 20 CRs of the BOX2560 (solid lines) and of the BOX1280 (dashed lines). Right panel: The 
mean (thick lines) and the mean plus/minus one standard deviation of the distribution of CRs (BOX2560 solid lines, BOX1280 dashed lines). 



WF (MB corr.) 


WF(no 

MB corr.) 

R 

Fbulk 

1 

b 

Vbulk 

1 

b 

10 

536 

282 

27 

513 

284 

27 

50 

250 

290 

18 

249 

295 

18 

100 

270 

283 

19 

274 

287 

19 

150 

239 

284 

11 

239 

288 

13 


Table Cl. The amplitude and direction (in Galactic (/, b) coordinates) of the bulk velocity for the WF velocity field calculated for Malmquist bias corrected 
data (denoted by ’MB corr’.) and for the raw uncorrected data (denoted by ’no MB corr.’). (Amplitude is given in km and distances in h“^Mpc.) 
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Figure Dl. Comparison of the mean (WF) and scatter (CRs) around the mean of the bulk velocity constructed from the CF2 (solid lines, black) and the CFl 
(dashed lines, blue) datasets. The case of unconstrained realizations (dotted lines, red) as well. The random realizations used to calculate the mean and cosmic 
variance are the same ones used to construct the CFl and CF2 constrained realizations. 
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